Robustness of Binary Black Hole Mergers in the Presence of Spurious 

Radiation 
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We present an investigation into how sensitive the last orbits and merger ol binary black hole 
systems are to the presence of spurious radiation in the initial data. Our numerical experiments 
consist of a binary black hole system starting the last couple of orbits before merger with additional 
spurious radiation centered at the origin and fixed initial angular momentum. As the energy in 
the added spurious radiation increases, the binary is invariably hardened for the cases we tested, 
i.e. the merger of the two black holes is hastened. The change in merger time becomes significant 
when the additional energy provided by the spurious radiation increases the Arnowitt-Deser-Misner 
(ADM) mass of the spacetime by about 1%. While the final masses of the black holes increase 
due to partial absorption of the radiation, the final spins remain constant to within our numerical 
accuracy. We conjecture that the spurious radiation is primarily increasing the eccentricity of the 
orbit and secondarily increasing the mass of the black holes while propagating out to infinity. 

PACS numbers: 04.25.ds,04.30.Db,04.70.Bw 



I. INTRODUCTION 



The coalescence of two black holes, long thought 
of as the "holy grail" of numerical relativity (NR) , 
is well on its way to being a solved problem. Many 
groups in NR have now demonstrated the ability 
to follow two black holes through several orbits 
[l| and their final orbits and merger to a single 
black hole 0, i, i, i, H, 0, S ■ From the first pub- 
lished waveform of equal-mass, non-spinning bi- 
nary black hole (BBH) coalescence, the simplicity 
of the waveform's dependence on time has been 
noted. Comparisons amongst the groups in NR 
have demonstrated a remarkable agreement to the 
solution of the BBH problem. A common aspect 
in all numerical relativity BBH evolutions is the 
presence of spurious radiation in the initial data. 
In this paper, we present a study on how the stan- 
dard equal-mass, quasi-circular BBH system re- 
sponds to the presence of spurious radiation that 
has been added in a controlled manner and map 
that response as a function of the radiation's initial 
conditions. Our intent is to determine how much 
junk radiation the system can handle and how the 
waveforms and the physical properties of the final 
black hole deviate from the standard BBH result. 



Several papers have compared BBH waveforms. 
One of the first papers to internally compare 
waveforms also demonstrated the first evidence of 
"universality" ^Qj in an equal-mass, non-spinning 
initial configuration. The paper demonstrated 
that differences in initial data characterized by 
a change in the initial orbital separation mani- 
fested as a time shift in the amplitude and phase 
of the gravitational waveforms. Once time-shifted, 
the waveforms were within 1% agreement over the 
merger and ringdown in |ro4'4|. We investigate 
the effect that the additional spurious radiation 
we add to the binary will have on this universal- 
ity. 

The first comparison of NR waveforms between 
several groups [10] includes the most popular 
methods used in the community to evolve BBHs 
covering excision with a hyperbolic formulation 
and moving punctures with the Baumgarte- 
Shapiro-Shibata-Nakamura (BSSN) formulation 
of the Einstein equation [1, |3| . The waveforms 
were in remarkable agreement once time-shifted, 
the largest differences, occurring at the beginning 
of the wave, being due to the spurious radiation in 
the initial data. A second, independent compari- 
son of waveforms from different methods was con- 
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ducted by Sperhake llj in which he compared a 
Kerr-Schild/excision evolution to a puncture evo- 
lution within the same code. An interesting issue 
to investigate is to what extent the spurious radi- 
ation in the initial data could cause differences in 
the merging time and thus affect waveform com- 
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parisons based on time-shifts to align the ampli- 
tude of the waveform. 

Most groups remove the initial burst from the 
waveform during post-processing of the data [HI, 
[Tsj . From the evolutions published, it appears 
that the spurious radiation that is present in the 
initial data is flushed out of the system within a 
crossing time, leaving the binary dynamics mostly 
untouched. There is still some concern about the 
impact that choices made in setting up the initial 
data for the evolutions, choices like conformal flat- 
ness, have on the waveforms. Studies have looked 
at different ways of choosing the freely specifiable 
part of initial data ^14] which reduce the ampli- 
tude of the spurious radiation, but these have not 
been extensively implemented in evolutions. 

In this paper, we test the robustness of the 
binary to effects of spurious radiation. We cre- 
ate a BBH system containing additional radiation 
with tunable initial energy initialized at the bi- 
nary's center of mass. We then evolve a series 
of equal-mass, non-spinning, quasi-circular BBH 
plus radiation spacetimes using the PSU numer- 
ical code that implements the Moving Puncture 
Recipe (MPR) ji, J, see [l^ for details about the 
PSU code. The initial data used to construct the 
BBH plus radiation spacetime is presented in Sec- 
tion IHII the results in Section IIVI and the conclu- 
sions in Section |Vl Our main result is that the 
presence of spurious radiation causes a hastening 
of the merger, thus plausibly accounting for the 
differences in merger times seen in the NR wave- 
form comparisons. Before describing our numeri- 
cal experiment, we present a back-of-the envelope 
calculation to build our intuition about this prob- 
lem. 



II. A NEWTONIAN PERSPECTIVE 

For illustrative purposes, we investigate the ef- 
fect a central pulse of energy might have on a bi- 
nary by studying a two-body orbit in Newtonian 
gravity with a stationary mass placed at the or- 
bit's center of motion while the bodies are at their 
apocenter. The addition of the third mass at the 
center of the Newtonian binary affects the orbit by 
deepening the potential in which the binary sits. 
We solve the problem using the standard central 
force solution to the two-body problem with the 
new potential. We assume r = initially since this 
is set in the parameters to the initial data solver. 
Letting m be the masses of the black holes and 



rrijju be the equivalent mass of the third body, we 
write the ratio of the final eccentricity, e', to the 
original eccentricity as 



e' (l + 2/)2 



^^^(4-2,V,d)/ + 4/2^ (1) 



where / = mw/m is the fractional mass, j ~ l/fi 
is the angular momentum per unit reduced mass, 
and d is the initial separation of the binary. 

This simple calculation indicates that the eccen- 
tricity increases for sufficiently small eccentricities. 
For the binary parameters studied in this paper, 
the eccentricity invariably increases for e < 0.88. 
Although the black holes in our BBH evolutions 
are not far enough apart to allow a valid determi- 
nation of eccentricity, the trajectories are quasi- 
circular enough for the eccentricity to be low. This 
illustrates that we can expect the addition of radi- 
ation into the studied system to cause the binary's 
orbit to become more elliptic. 



III. INJECTING RADIATION INTO A 
BBH EVOLUTION 

We inject radiation into the standard, equal- 
mass, non-spinning, quasi-circular BBH evolution 
during the setup of the initial data. The initial 
data for the evolution is constructed via the punc- 
ture method [l^ using the single-domain spectral 
method code developed by Ansorg et al. [I4I that 
uses a conformally flat prescription to solve the 
constraints. We have two building blocks for the 
data: 1) the quasi-circular BBH and 2) the tun- 
able radiation. The BBH data is set-up using the 
input conditions for the Baker et al. [!|] Rl run of 
two equal-mass irrotational black holes in quasi- 
circular orbits such that the metric is conformally 
flat. The details of the Rl initial data are given in 
the first row of Table [T] and a convergence study 
was done in llal. 



A. The Teukolsky-Nakamura Wave 

The tunable radiation is given by an even par- 
ity, quadrupolar gravitational wave: the linearized 
solution to a perturbation on Minkowski space- 
time expanded over the modes of the Matthews 
tensor spherical harmonics [l^. This radiation 
was first derived by Teukolsky and is typi- 
cally known as a Teukolsky wave. This tensor is 
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simply the general traceless-transverse solution to 
the linearized Einstein equation. Nakamura (20| 
later wrote out the solution for general i and m. 
Instead of using that tensor as a metric pertur- 
bation, he used it as an extrinsic curvature per- 
turbation. We implement the Nakamura version 



Note that the TN solution lets us choose the 
radial dependence in the form of ingoing and out- 
going functions which we have chosen to be the 
same symmetric functional form, F(u). F(u) is 
given by an Eppley packet [2l[ 

F{u) = Awe-"'/'"' , (5) 

where u ~ t ± r . The Eppley packet is localized 
and smooth with the factor of u present so the 
wave is more regular at the origin. Later in the 
paper we investigate a cosine modulation of this 
function. We can choose the location, mode con- 
tent, strength, and radial dependence of the in- 
jected radiation. Since the extrinsic curvature is 
real and the spherical harmonics are complex, we 
take only the real part of the Xim tensor result- 
ing in a superposition of m and — m modes in the 
TNW. 



of the Teukolsky wave, herein called Teukolsky- 
Nakamura waves (TNWs), in order to satisfy the 
condition of a conformally flat metric imposed by 
the puncture method. 

The Teukolsky-Nakamura (TN) extrinsic curva- 
ture tensor is given by 



(2) 

(3a) 
(3b) 

(3c) 
(3d) 

(4a) 
(4b) 

I 

B. BBH+TNW 

Our initial data, the spatial metric and extrinsic 
curvature, are given by 

mj = i^'^Vtj, (6a) 
K., = V-^(i^^+i^.^), (6b) 

where rjij is the flat spatial metric and ijj is the so- 
lution to the Hamiltonian constraint under York's 
conformal approach [2^ . The extrinsic curvature 
Afj^ is the Bowen and York solution [2^ to the 

momentum constraint and Aj,'^ is the TN tensor. 
Notice that because the momentum constraint is 
linear in the extrinsic curvature, the superposition 
of the extrinsic curvature is also a solution of the 
momentum constraint. As a test we evolved the 
above initial data with a vanishing Bowen- York 



(aimYi 
gimYlm + flmWlm flmXlm 
{gimYlm - flmWlm) siu^ I 

where the coefhcients a/,„, him, fim, and gim are functions only of r and t as follows 
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tensor, i.e. Minkowski background. For the i — 2, 
m = ±2 case, we found that the Arnowitt-Deser- 
Misner (ADM) angular momentum calculated on 
the initial spacetime is zero to within machine er- 
ror. The angular momentum of the BBH+TNW 
is therefore independent of the TNW to our nu- 
merical accuracy. 

C. Configurations 

The simplest geometry to add additional spuri- 
ous radiation to our BBH initial data is to center 
a wave pulse at the origin. We typically choose 
an ^ = 2, TO = 2 mode, the dominant mode for 
gravitational radiation from a BBH system, and 
vary the amplitude, A, and width, cr, of the Epp- 
ley packet. The values oia/M are chosen from the 
set {0,3,4,5,6} and those of the amplitude A/M^ 
from the set {0,0.1, 0.5, 1.0, 1.5} where the dry 
Rl BBH spacetime is recovered when A/M^ — . 
Fig. [T] shows the shape of the wave in one of the 
components of the traceless-transverse extrinsic 
curvature, Aij = ^~'^Aij, along the coordinate 
axis intersecting the two black holes, modulated 
by the inverse square of the conformal factor, for 
A = IM^ and cr = 3M, 4M, 5M, 6M. 
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FIG. 1: Comparison of Axx initially between various 
Eppley packet widths, a, for a wave amplitude of IM'^ . 
This is from the binary black hole initial data where 
Axx vanishes at the punctures. 

When adding the TNW to the spacetime we 
wanted to keep the initial black holes unaltered. 
We chose to keep the apparent horizon (AH) 
masses constant independent of the additional 
wave content. In practice the AH masses varied 
by as much as 0.04% from the dry Rl run due 
to insufhcient parameter accuracy. The momenta 
remain constant as parameters to the initial data 



solver, and the ADM angular momentum differs 
by at most 0.001% from the dry Rl run. The sec- 
ond column of Table |l] lists the ADM energy of the 
spacetime for each wave choice. We note that the 
waves increase the ADM energy from a negligible 
10^^^% to a significant 8.9%, which scales empir- 
ically as -Eadm — A^ /a^ . The proper separation 
between the black holes changes from the dry case 
oi L = 9.94 to a maximum of i = 10.23. In the 
most extreme case, the wave having A — 1.5M'^ 
and CT = 3M, we have added almost 9% additional 
energy into the BBH system and effectively moved 
the black holes apart by 0.29M. The impact of 
these differences in initial data on the binary evo- 
lution are discussed further in the next section. 

Pumping energy into the system while holding 
the coordinate separation and angular momentum 
constant necessarily means that we are changing 
the binding energy of the system. To study this 
change we map the effective potentials for each 
BBH-I-TNW case. We do this by repeatedly solv- 
ing the initial data with incremented separations 
while holding the individual AH masses and to- 
tal ADM angular momentum fixed. We calcu- 
lated the quantity Ei, = Eadm — "T'Ah.i — m'AH,2 
for each spacetime. The wave itself adds to the 
ADM energy and must be subtracted in this cal- 
culation; but, as we're only interested in the rel- 
ative shapes, we can look at the relative binding 
energy. Eh — -Bh^min- For the waves with a cr of AM 
the binding energy per unit reduced mass is plot- 
ted in Fig. [2] with a vertical line indicating the ini- 
tial coordinate separation of the black holes. We 
can immediately see that the dry "quasi-circular" 
Rl case has some non-zero eccentricity as the im- 
posed separation does not lie at the minimum of 
the curve. We also observe a shift of the minimum 
inwards as the wave strength increases. Since the 
coordinate separation in the parameter search is 
held fixed for the evolved initial data, the location 
of the system along the binding energy curve with 
respect to its minimum is sufficient to see that the 
eccentricity of the orbit is likely increasing. Un- 
fortunately, the separation is too small to get a 
reliable measure of the eccentricity. 



IV. EVOLUTIONS OF THE BBH+TNW 

Our simulations of the BBH+TNW initial data 
are summarized in Table HI where the first row 
corresponds to the BBH system without added ra- 
diation, the dry Rl run. This will be our control 
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FIG. 2: We plot the effect of the wave on the binding 
energy per unit reduced mass in the initial data. The 
potentials were calculated by solving the initial data 
using Ansorg's code for various separations while keep- 
ing the individual AH masses and total ADM angular 
momentum fixed. 



case. We systematically evolved each BBH+TNW 
spacetime, varying A and a of the TNW. We di- 
vide the results from evolving these simulations 
into four subsections: the main result concerning 
the merger time in § IIV Al the dynamical and ra- 
diated quantities from our runs are in § IIVBI the 
final spacetime quantities in ? IIV CI and a compar- 
ison of the gravitational waves by time-shifting in 
§[IVDl 



A. Merger Time 

The main result of adding gravitational radia- 
tion to our BBH evolution is to hasten the merger 
of the black holes. With increasing Eadm, the bi- 
nary invariably merges faster. The sixth column 
in Table [T] lists the differences in merger times be- 
tween the dry Rl and the BBH-I-TNW runs given 
by AT ^ (To - To^dry)/M . The time, To, is given 
in units of the total, initial AH masses of the black 
holes and evaluated at the peak amplitude of each 
waveform extracted at a radius of 75M . The use 
of the waveform peak variation as a measure of the 
change in merger time agrees within a few percent 
in AT to the variation in the time it takes for the 
punctures to be separated by one grid spacing. 

Figs. [3] and S] show the change in merger times 
from the perspective of constant wave amplitude 
and constant pulse width. We can see that there is 



A/M'-' a/M 


Badm Jj:ad_ AT/M M, j. 


0.0 
0.1 3 
0.5 3 
1.0 3 
1.5 3 


0.9957 0.0359 0.273 0.0 0.9599 0.682 
0.996 0.0363 0.273 -0.6 0.9600 0.683 
1.007 0.0451 0.271 -16.4 0.9609 0.682 
1.037 0.0708 0.263 -56.1 0.9635 0.686 
1.084 0.1058 0.244 -88.4 0.9696 0.693 


0.1 4 
0.5 4 
1.0 4 
1.5 4 


0.996 0.0360 0.274 -1-1.5 0.9596 0.682 
0.999 0.0385 0.272 -4.7 0.9604 0.682 
1.007 0.0463 0.272 -14.3 0.9603 0.683 
1.021 0.0589 0.270 -31.8 0.9607 0.683 


0.1 5 
0.5 5 
1.0 5 
1.5 5 


0.996 0.0360 0.273 -fO.2 0.9599 0.682 
0.997 0.0369 0.273 -0.1 0.9599 0.682 
1.000 0.0399 0.272 -4.6 0.9603 0.683 
1.005 0.0448 0.272 -7.8 0.9599 0.686 


0.1 6 
0.5 6 
1.0 6 
1.5 6 


0.996 0.0359 0.273 -fO.2 0.9599 0.682 
0.996 0.0364 0.273 -fO.7 0.9599 0.682 
0.998 0.0377 0.273 -0.3 0.9601 0.682 
1.000 0.0399 0.272 -2.2 0.9602 0.682 



TABLE I: The first two columns are the parameters of 
the TNWs followed by the ADM energy of the initial 
spacetimes. Column 4 and 5 give the fraction of the 
ADM energy and angular momentum radiated over 
the simulation. Column 6 is the change in merger time 
calculated by the shift in extracted waveform peak in 
units of the total Alf mass in the initial spacetime. 
Column 7 lists the final mass and Column 8 the final 
spin, ii — af/Mf of the black hole. 



a strong dependence on the width of the pulse as 
well as the amplitude. Some cases show a positive 
value for AT; however, these are all equal to zero 
within the errors. For all the A = IM^ waves that 
have non-zero merger time we found an approxi- 
mate power law relation between the width of the 
pulse and the change in merger time: 



AT(yl = IM^) cx cr" 



(7) 



A more general look at the change in merger times 
is found in Fig. [5] Given our estimated error bars, 
significant changes in merger time occur when the 
TNW has increased the initial ADM energy of the 
spacetime by about 1% compared to that of the 
dry Rl. 

To isolate how much of AT is due to the ad- 
ditional spurious radiation introduced and how 
much is due to other factors, we perform a series 
of tests. We focus on the most significant sources 
of errors, namely the resolution of our grid, wave 
extraction radius, the change in proper distance 
in setting up the initial data, and the change in 
mass of the two black holes. We will look at each 
of these factors and assess their individual contri- 
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FIG. 3: Changes in merger times compared to the dry 
Rl run as a function of packet width with estimated 
error bars. 



FIG. 5; Changes in merger times compared to the frac- 
tional change in initial ADM energy of the spacetime 
with estimated error bars. 
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FIG. 4: Changes in merger times compared to the dry 
Rl run as a function of wave amplitude with estimated 
error bars. 



bution to AT. 

1. Resolution: The finest resolution for the sim- 
ulations we present in Table |T] is M/38.4. We 
check the error due to the resolution by repeat- 
ing several cases with finest resolutions of M/44.8 
and M/51.2. We ran convergence tests on the 
strongest wave {A — 1.5M^, a — 3M) and the 
weakest wave (A = O.IA/''', a = 6Af). A third, 
medium, wave with A = 0.5M3 and a = AM 
was run at just one more resolution, Af/44.8. For 
all three cases, the merger time decreased. The 
merger time of the compact wave decreased lit- 
tle for a total of 0AM over the three resolutions 



while the diffuse wave decreased more drastically 
for a total of about 0.5M. The medium case had 
a 0.2A/ difference between the two resolutions. 

2. Extraction Radius: The next source of er- 
ror is wave extraction radius. In NR, waveforms 
are usually calculated in terms of the Newman- 
Penrose scalar, \I'4(t, a;, y, z), which are extracted 
on a sphere at a finite radius some distance 
from the source, then expanded into angular 
modes via the spin- weighted spherical harmonics, 
-2Yem{8,(f>)- With a proper choice of tetrad, this 
scalar is a measure of outgoing gravitational ra- 
diation. There has been recent work investigating 
the effects the choice of extraction radius can have 
on the correctness of the waveform [2^. As 
the extraction radius increases, the errors caused 
by an incorrect tetrad and finite distance dimin- 
ish. While it is still an open question whether 
or not there are observable effects from the meth- 
ods groups currently use to extract the waveforms, 
the methodology of the extraction is not thought 
to contaminate the waveform. An indication the 
appropriate tetrad is being approached is that 
the waveform amplitude scales as l/r, which we 
tested. To get a rough estimate of the errors due 
to extracting at a finite radius, we compute AT 
using radiation extracted at 30M and extracted 
at 75M. In Fig. [6] we plot the amplitude of the 
dominant waveform mode, |^'4'^|, extracted at the 
two radii for both the dry Rl run and one where 
the merger time changed significantly. The merger 
time shift changed by no more than 0AM between 
the two extraction radii. 
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FIG. 6: Sample comparison of waveforms extracted at 
different radii. Plotted are the waveforms for the dry 
Rl run and A = 1M^ a = 5A/ run extracted at 30M 
and 75M. 



3. Black Hole Mass: Aside from unphysical 
sources of error, the small differences in the ini- 
tial data also change the merger time. Though we 
kept the initial AH masses nearly identical, there 
is still a variation of up to 0.04% compared to the 
dry Rl run. While conducting the research for this 
paper, we found that a 0.14% change in initial AH 
masses of the punctures resulted in a change in 
merger time of 6.7M. Although we do not include 
simulations with such a large deviation of masses, 
we used this knowledge and assumed the change 
in merger time was linear in the change in initial 
mass to estimate an error. 

4-. Proper Separation: As mentioned in § IIIIl 
the presence of the additional gravitational radia- 
tion also increased the proper separation, L, from 
the dry Rl 9.94M by up to 2.9%. We studied the 
effect of this change by evolving two BBH space- 
times with the same initial masses and angular 
momentum but increasing the coordinate to yield 
proper separations of 10.24M and lO.lOM. The 
merger time changed by at most 2.1M. We as- 
sumed a linear relationship between the AT and 
AL in estimating the errors from this source at 
each data point. 

The error bars presented in our figures are cal- 
culated by adding all the errors in quadrature: 



(8a) 



AL Am 



T, and is the error in T for the dry Rl run. 
Srcs = 0.75 and Etet = 0.4 are the largest mea- 
surements for the resolution and extraction radius 
errors. The accumulated errors do not account for 
the observed AT when i?ADM/^^Ri > 0.01 and 
we note that the errors grow as the amplitude in- 
creases and the width decreases, most notably the 
errors associated with changes in the irreducible 
masses. 

The parameter space of adding spurious radia- 
tion is large. In Table |TT1 we present the results 
from a few evolutions outside of our main param- 
eter survey. The junk radiation present in the ini- 
tial data of a typical BBH simulation may not be 
well represented by an £ = m = 2 mode. Similarly 
the effect of the junk radiation might be sensitive 
to the wavelength of the pulse. In order to test how 
important a modulation in the frequency might be 
to our conclusions, we briefly investigated an Ep- 
pley packet modulated with a cosine wave, given 

by 



F{u) = Acos(fcM)we^"'/'"' 



(9) 



where Eat is the error in AT, Et is the error in 



This modulation adds an extra parameter control- 
ling the wavelength of the perturbation. We ad- 
justed the amplitude of the wave to keep the en- 
ergy approximately comparable to our standard 
runs. The resulting simulation merger time dif- 
fered from the unmodulated packet by less than 
IM in T, well within error bars. While this is 
still an avenue open to investigation, we concluded 
that the modulation was not affecting the results 
enough to warrant an additional parameter in our 
survey. We also conducted a test of the geometry 
of the wave by initiating a pulse with an £ = 2, 
TO = mode. Again we changed the amplitude 
so that the energy in the wave was approximately 
constant and found that there was a change in 
merger time of IM compared to the £ = 2, m = 2 
simulation, again within error bars. This points 
towards the wavelength and angular dependence 
of the pulse being secondary to the additional en- 
ergy in determining the effect of the pulse on the 
merger time. 

Finally, to make a stronger connection to the 
junk radiation being associated with each punc- 
ture, we added two identical waves centered at 
each of the black holes rather than at the center. 
Compared to the same wave initiated at the cen- 
ter, the dual waves added almost twice the energy 
and almost doubled the change in merger time, 
which is consistent with the center-of-mass TNW. 
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TABLE 11: Overview of the odd runs. The left four 
columns are the wave parameters, followed by the frac- 
tion of ADM energy radiated and the fraction of the 
ADM angular momentum radiated. The final masses 
are given in the 7th column followed by the merger 
time change as derived by the peak of the waveform 
extracted at 75M. 



B. Dynamical and Radiated Quantities 

We now investigate the effects of the TNWs 
on the radiated quantities derived from the wave- 
forms. We calculate these quantities from the 
Weyl scalar ^4 assuming, as in the calculation 
of the waveform, the fiducial tetrad of Baker et 
al. [26j . A summary of the quantities obtained 
from ^4 are listed in Table HI The fraction of 
the initial ADM energy radiated was calculated 
across a detector at 40 Af. As expected, the ra- 
diated energy increases with the strength of the 
wave. When A < 0.5AI^ and a > 5M, there is 
no measurable difference between the BBH-I-TNW 
and the BBH cases within numerical error. For 
those cases, we can only conclude that the energy 
in the wave propagates out without a measurable 
interaction with the black holes. Similarly a trend 
emerges as we increase A for each cr, which corre- 
sponds to increasing Eadm ■ The radiated angular 
momentum consistently decreases. 

To look at the interaction of the TNW with the 
black holes as it propagates out, we study the radi- 
ated energy and angular momentum as functions 
of time. In Fig. [7] we plot the energy radiated 
across a detector at radius r = 40M. We see the 
energy grows from a time of 40M to around 80M 
as the initial burst of spurious radiation passes the 
detector. After this burst of energy the remaining 
energy radiated is approximately 0.035£^adm,ri 
and is almost uniform across the various cases. 
From this we can see that most of the energy in- 
troduced in the spacetime is quickly flushed out 
of the system, leaving a system which radiates a 
further amount of energy that is independent of 
the junk radiation. 

All the evolutions started from spacetimes with 



equal Jadm since the TNW does not add angular 
momentum to the BBH spacetime. In column five 
of Tabic [J and in Fig. [51 we see that the amount 
of angular momentum radiated across a detector 
located at 40M is independent of the wave with 
some numerical error. The difference in Jrad be- 
tween the runs lies in when the system radiates 
the angular momentum. This is better seen in 
Fig. [9] where we present the angular momentum 
flux across the sphere at r = 40M and in Fig. \10\ 
a close-up of the initial part of the data. The spu- 
rious radiation is transporting extra angular mo- 
mentum as it is flushed out. 

Figs. [7] and [5] show some of the TNW propagat- 
ing out at the same time as the spurious radiation 
is flushing out. One of the consequences of the 
wave traveling away from the center of the orbit 
is that it can interact with the black holes and 
potentially increase the mass of each black hole 
during the early part of the inspiral. Table IIIII 
documents how the AH masses change as a func- 
tion of A and a for the stronger TNW cases. The 
black-hole mass is calculated using an AH tracker 
[27I I. The Am is a measure of the change in the 
mass of each black hole up to 50M, such that 
Am = m{t — 50M) — m(< — 0) . The change in 
the initial ADM energy compared to the Rl run 
is given by AM = E'adm - -Eadm.ri at t = . 
We use this estimate of the differences in the mass 
of the spacetime between Rl and the rest of the 
runs to compute a naive estimate of the total frac- 
tion of energy absorbed by both black holes. Up 
to 7.89% of the extra ADM energy is observed to 
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FIG. 7: Energy radiated across a sphere of radius r = 
4:0M as calculated from the Weyl scalar ^'4. 
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FIG. 8; Angular momentum radiated across a sphere 
of radius r = 40M for a of 3AI as calculated from the 
Weyl scalar *1'4. 



be absorbed by the black holes during the first 
50M of the simulation, 3.9% by each black hole. 
The actual amount absorbed depends strongly on 
the width of the wave: the narrow, strong pulses 
are more readily absorbed than the weak, diffuse 
pulses that extend beyond the black holes in the 
initial data. In the weaker cases the change is 
barely visible above the noise in the AH mass cal- 
culation, in the stronger cases it is unmistakable. 

To assess how important the absorption of en- 
ergy by the black holes during the evolution is to 
the changes in merger time and radiated angular 
momentum, we refer to our discussion about the 
sensitivity of the merger time to a change in the 
initial AH masses in § IIV Al Given a change in 
mass of the individual black holes of 0.14%, the 
merger time changed by 6.7M. In setting up the 
initial data, we do not allow the AH masses to 
change more than 0.04%. The amount of absorp- 
tion measured during the evolution is as much as 
3.9%; and, therefore, the increase in mass may be 
accounting for some, although not all, of the ef- 
fects of the TNW. The outliers, the cases of most 
extreme merger times, merge so quickly that dif- 
ferentiating the burst of spurious radiation and 
region of pure inspiral is difficult. We also tried 
normalizing the time axis by the total AH masses 
after the wave has passed rather than that from 
the initial spacetime. This changed AT by no 
more than IM so the choice of normalization does 
not account for the observed difference in merger 
times. 



FIG. 9: Flux of angular momentum across a sphere of 
radius r = 40M for a of 4M as calculated from the 
Weyl scalar 'I'4. 
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FIG. 10: Closeup of angular momentum flux across a 
sphere of radius r = 40Af concentrating on what the 
spurious radiation carries. 



C. Final Spacetime 

One of the important products of a BBH coa- 
lescence to relativists and astrophysicists are the 
final black hole's mass and spin. The final black- 
hole masses and spins are presented in the last two 
columns of Table [J To compute the final mass, 
Mf, we use energy conservation arguments by cal- 
culating the difference between the ADM energy 
and the radiated energy as calculated from the 
Weyl scalar ^'4 . The final spin, jj = aj/Mf, is 
calculated by finding the complex ringdown fre- 
quency in the £ = 2, m = 2 mode and using the 
numerical Kerr frequencies given in Table II of Ap- 
pendix D in Berti et al [28|] to find the correspond- 
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TABLE III: Change in AH mass compared to the dif- 
ference in initial ADM energy for the stronger waves. 
Am is the change in a single black hole AH mass over 
the first 50Af , AAf is the additional ADM energy com- 
pared to the dry Rl run. Column 5 is the fraction of 
the extra ADM energy absorbed by both black holes 
combined, and the last column is the ratio of the final 
black hole mass to the total AH mass after the wave 
has been absorbed. 



ing spin parameter. This method agrees within 
stated errors to inverting the fit of Eq. (E2) of the 
same paper. Given the strong dependence of the 
spin on the damping time, we limit ourselves to 
the real part of the complex frequency and com- 
pare this to a separate spin calculation using the 
isolated horizon framework [29] where possible. 

From the values of Mf and jj listed in Table [H 
we can see that the final spins are constant within 
numerical accuracy and the final masses do not 
vary strongly with A and a. The trend is an in- 
crease in the final mass with increasing £^adm be- 
coming noticeably greater than our numerical er- 
rors for the four largest cases, Mf > 0.963 . From 
this we can see that the narrower pulses not only 
have more energy, but they also interact more ef- 
ficiently with the black holes. Being more readily 
absorbed by the punctures, they increase the in- 
dividual masses and thus the final mass. The last 
column of Table IIIII shows the ratio of the final 
mass to the total AH mass once the wave has in- 
teracted with the inspiraling black holes. We see 
the ratio is roughly constant, implying that ap- 
proximately 5% of the initial AH mass is radiated 
away if we include the wave energy absorbed by 
the black holes. The exception is the most ex- 
treme wave where the black holes merged before 
all the spurious radiation has been absorbed into 
the AH. This would underestimate the AH growth 
and thus overestimate the value of the ratio. The 
change in final mass agrees within numerical error 
to the change in total AH mass after the wave has 



D. Alignment of Amplitude and Phase 



As stated in the introduction, a common 
method to compare waveforms is via a time-shift 
of the amplitude of each waveforms such that their 
peaks overlap, the result of which is shown in 
Fig. [TT] We can see that the waveforms overlap 
very well after the merger. The only noticeable 
difference is for the strongest TNW we evolved 
on the BBH system, the A = l.SM^, a = 3M 
case, where we find the largest difference in the 
final black hole compared to the dry Rl run. We 
can also see residual contamination of the merger 
portion of the waveform by the spurious radiation 
due to the binary merging so quickly. Similarly, 
we shift the waveform phase such that they over- 
lap at T = Tpoak in Fig. [T^ The agreement in the 
phase's slope during ringdown is further confirma- 
tion that the mass and spin of the final black hole 
are not significantly altered. 

The waveform overlap in the merger regime con- 
tinues to start before the merger, as seen in jl(| . 
as long as the spurious radiation does not con- 
taminate this region of the waveform. This al- 
ludes to the relatively simple form of the merger 
waveform seen in all the various situations cur- 
rently tested. As long as the spurious radiation 
is not strong enough to noticeably alter the final 
black hole, the merger portion of the waveform re- 
mains essentially unaltered and the contamination 
to the system predominately results in the change 
in merger time and thus a time-shift of the wave- 
form. 



V. CONCLUSIONS 

In this study, we simulated an equal-mass, non- 
spinning BBH system through its last orbits, 
merger and ringdown. The system is perturbed 
by the systematic addition of spurious radiation 
in the form of a Teukolsky-Nakamura gravitational 
wave at the binary's center of mass. The initial en- 
ergy of the wave is tunable, specified by the ampli- 
tude and width of the radiation; in addition, the 
initial angular momentum was fixed for the entire 



11 





0.08 




0.07 


T3 

<D 


0.06 




0.05 








0.04 








0.03 


O 

l-H 


0.02 




0.01 









Rl 

Afi.5M-i,a=4M 
1.5 ]V^,a=3M 



-40 



-20 





(T-To)/M 



20 



40 



FIG. 11: Comparison of for all runs shifted such 

that the peak amplitudes before ringdown coincide. 
This lets us compare the relative damping times of 
the ringdown and thus the properties of the final black 
hole. We also note the agreement for about 50M be- 
fore ringdown as well. Though the legend only labels 
the distinguishable cases, all the runs are contained in 
the figure. 
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FIG. 12: Comparison of arg(*4'^) for all runs shifted 
such that the peak amplitudes before ringdown coin- 
cide and shifted vertically. This lets us compare the 
relative frequencies of ringdown. Though the legend 
only labels the distinguishable cases, all the runs are 
contained in the figure. 



sequence of runs. The binaries that contain the 
extra radiation invariably merge faster than those 
with no additional radiation. 

In addition to the main result of decreasing 
merger time, some changes to the radiated quan- 
tities and the final black hole were measurably 



above numerical error. This occurred once the ad- 
ditional energy provided by the TNW was equal 
to or greater than 1% of the dry BBH spacetime. 
As the TNW propagated out of the center, ap- 
proximately 4% of additional ADM energy was 
absorbed by each black hole. In that strong-wave 
case, it was not possible to make an accurate mea- 
surement of the mass of the enlarged black holes 
before the plunge of the binary. The final spins of 
the black holes, however, remained unaffected by 
the gravitational radiation for all but the strongest 
case (A = 1.5AP, a = 3M). The constant black- 
hole spin is consistent with the wave slightly in- 
creasing the eccentricity of the orbit for small ec- 
centricities [13, HH . We also observed a decrease 
in the radiated angular momentum with increas- 
ing TNW strength. 



We conjecture, based on the change in the initial 
binding energy of the BBH+TNW systems and 
Newtonian back-of-the-envelope calculation, that 
the spurious radiation increases the eccentricity of 
the original orbit. Unfortunately, the separation of 
the black holes was not large enough to enable a 
reliable calculation of the eccentricity. The merger 
time is very sensitive to the increase in individual 
black-hole masses via wave absorption; however, 
this was not enough to account for the observed 
change in the time of merger even when ignor- 
ing the strongest wave case. The combined effects 
of increasing the individual black-hole masses and 
the eccentricity of the orbit caused the binaries to 
merge faster with increasing energy. 



One of the conjectures in the literature is that 
the spurious radiation, intrinsic to the construc- 
tion of initial data for BBH evolutions, is flushed 
out of the simulation within a crossing-time and 
does not effect the radiation or the binary. We 
can relate the results of this study to other BBH 
evolutions by looking at the early changes in AH 
mass as well as how much energy leaves the sys- 
tem in the burst of spurious radiation. For the dry 
Rl run, the energy radiated in the initial pulse is 
9 X 10~^i?ADM,Ri- We find that there is negligible 
effect on the merger time at that level. Our re- 
sults indicate that the spurious radiation present 
in initial data sets is unlikely to cause dramatic de- 
partures from the true BBH solution and therefore 
we can state that the simulated merger is robust 
to the presence of spurious radiation. 
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